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ABSTRACT 

We apply a new non-parametric technique to reconstruct, with uncertainties, the projected 
mass distribution of the inner region of Abell 2218, using combined strong and weak lensing 
constraints from multiple-image systems and arclets with known redshifts. The reconstructed 
mass map broadly resembles previous, less detailed, parametric models, but when examined in 
detail shows several sub-structures, not necessarily associated with light but strongly required 
by the lensing data. In particular, the highest mass peak is offset by ~ kpc from the main 

light peak, and projected mass-to-light in the directions of different cluster galaxies varies by 
at least a factor of 10. On comparing with mass estimates from models of the X-ray emitting 
gas, we find that the X-ray models under-predict the enclosed mass profile by, at least a factor 
of 2.5; the discrepancy gets worse if we assume that mass traces light to the extent allowed by 
the lensing constraints. 

Subject headings: Dark matter - galaxies: clusters: individual (Abell 2218) - gravitational 
lensing: weak & strong lensing - X-ray 



1. INTRODUCTION 

Galaxy clusters are thought to be the largest and 
most recently assembled gravitationally bound enti- 
ties in the universe and, thus, an accurate determi- 
nation of their masses is of utmost importance to 
full understanding of the formation and evolution of 
cosmic structures, mapping of dark matter in large 
scales and for constraining cosmological parameters. 
Traditional methods for inferring mass distributions 
in clusters are (i) dynamical methods, in which the 
observed line-of-sight velocity distribution of the lu- 
minous cluster galaxies is used in conjunction with 
the virial theorem, and (ii) X-ray methods, in which 
the investigation of the X-ray emission of the intra- 
cluster hot gas is used to trace the cluster potential. 
Both these traditional methods depend on restric- 
tive assumptions about the geometrical and dynam- 
ical state of clusters. A more sophisticated method, 
which proved to be very reliable and independent of 
any prior assumptions, is gravitational lensing. 

Theoretically, the above three methods have to 



yield the same cluster masses if the clusters are dy- 
namically relaxed, but recent literature shows a lot 
of discussion of and disagreement of cluster masses 
derived from observed velocity dispersions, observa- 
tions of X-ray emitting gas and gravitational lens- 
ing (Miralda-Escude and Babul 1995, Wu and Fang 
1997). Early studies based upon few selected clus- 
ters claimed a mass estimate discrepancy of at least 
~ 2 — 3 among them. This discrepancy can well be 
attributed to projection effects, non-thermal pressure 
or the specific cluster is still in the formation era. Ev- 
idently, the problem is relevant to the precise deter- 
mination of mass distribution and to the dynamical 
evolution of those clusters. For clusters which are still 
in an ongoing merging phase, i.e., not in hydrostatic 
equilibrium, the hot gas does not follow the gravi- 
tational potential of the cluster and thereby the X- 
ray cluster mass is uncertain and should be different 
from the gravitational lensing-derived mass and/or 
even the virial mass. 
Within gravitational lensing, there are two distinct 



* Based on observations made with the NASA/ESA Hubble Space Telescope, obtained from the data archive at the Space 
Telescope Science Institute. STScI is operated by the association of universities for research in Astronomy, Inc. under the NASA 
contract N-AS 5-26555. 
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methods for inferring masses from lensing observa- 
tions; (i) the parametric model-fitting method for 
the strong lensing region (giant arcs and multiple im- 
ages), and (ii) statistical distortion method for the 
weak lensing regions (weakly distorted single arclets) . 
The parametric model-fitting starts from the inner 
most regions of galaxy cluster, where lensing infor- 
mation is richest and extends outwards. It requiers 
fitting parametrized profiles to cluster-sized galaxies 
(Kneib et al. 96 [hereafter KESCS96]). It works ef- 
ficiently with one or two multiple image systems but 
becomes difficult to implement with several multiple 
image systems spanning outwards and hence the fit- 
ted model becomes non-unique. On the other hand, 
the statistical distortions method starts from regions 
where the distortions are very tiny and weak, i.e outer 
regions of galaxy clusters, and goes inwards. It re- 
quires averaging statistically the ellipticities of back- 
ground galaxies over patches of sky, typically ~ 20" 
in size, to trace the cluster mass distribution from its 
shear field (Kaiser and Squires 1993 [hereafter KS93], 
Kaiser 1995, Schneider and Seitz 1995). Methods, 
which are based on KS93 algorithm, suffer from a 
global invariance transformation known as the "mass- 
sheet-degenracy" and also cannot be used to probe 
the inner regions of rich clusters, because the lens- 
ing properties change very rapidly over that sampling 
scale. 

In intermediate regions, with no occurrence of mul- 
tiple images but having highly distorted arclets, nei- 
ther of the above methods can be used. However, 
several attempts to reonstruct the mass distribution 
in this region were made (Kaiser 95, Schneider and 
Seitz 1995, Seitz & Schneider 1996), but their for- 
mulation of the problem made it non-linear and still 
suffer from the so called mass-sheet degeneracy and 
boundary effects. 

Obviously, the problem of mapping the cluster 
mass distribution from gravitational lensing require 
an independent non-parametric method that can be 
used simultaneously in regions with varying lensing 
strength. The required method should disentangle 
the systematic effects that plagued the existing meth- 
ods encoded in forms of non-uniqueness. 

In this paper we present a new and general recon- 
struction technique that combines regions with vary- 
ing lensing strength the strong and overcomes the 
drawbacks of the existing methods. The technique 
is basically an extension of the non-parametric clus- 
ter inversion described by AbdelSalam ct al. (1998) 
[hereafter Paper I] for the strong lensing regime to 
incorporate the extra information encoded in the ob- 
servations of single/weakly distorted arclets. The ex- 
tension results in a method that combines strong and 
weak lensing data in a mass reconstruction and over- 



comes the mass sheet degeneracy. In general, the 
method is non-parametric and similar to the weak 
lensing (or statistical distortions) method of KS93, 
but our formulation of the problem is manifestly lin- 
ear in all the regimes and hence the problem becomes 
simpler and can thereby easily include the strong 
lensing regime. 

We apply the technique to the spectacular cluster- 
lens A2218 to reconstruct its projected mass distri- 
bution, with uncertainties, from combined strong and 
weak lensing. This cluster has previously been stud- 
ied through parametric lens modeling of strong lens- 
ing and arcs by Kneib et al. (1995) and KESCS96, 
non-parametric inversion from statistical distortions 
by Squires et al. (1996), and parametric modeling of 
statistical distortions by Smail et al. (1997). 

2. THE RECONSTRUCTION METHOD 

The method we will follow here is a free-form or 
non-parametric one that reconstructs a pixellated 
mass distribution. It is basically the same as in Pa- 
per I, except that (a) in this work we implement con- 
straints from weak as well as strong lensing whereas 
the earlier work implemented only strong lensing and 
described the extension to weak lensing, and (b) this 
paper uses Gaussian pixels whereas Paper I used 
square pixels — the difference is tiny for the results 
of either paper. Since Paper I already has a full de- 
scription, we will only summarize the technique here. 

2.1. Pixellization of the mass distribution 

The lens plane is divided into N x N pixels with 
inter-pixel distance a. The mnth pixel is a Gaussian 
tent with 'dispersion' a/2 and peak height n mn . That 
is to say, if the mn-th pixel is centered at mn it has 
a mass profile 
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Hereafter, we refer to a as the pixel size. We measure 
K mn in units of the critical density E C rit for sources 
at infinity. Thus the total mass is 
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The appropriately scaled arrival time of a light ray 
from a background source at an unlensed angular po- 
sition (3 via a point in the pixellated lens plane to 
the observer is 

T(0) = 1(0 - /3) 2 - ^ ]T Kmn^mniO), (2-3) 
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Thus, K m nipmn(6) represents the contribution of the 
mra-th pixel to the total gravitational potential of 
the lens. Its derivatives with respect to represent 
contributions to the bending angle and amplification 
components. 

See Appendix A for explicit expressions for the 
function ip mn and its derivatives. 

2.2. Lensing observables and constraint equations 

The basic approach of our technique is that the lens 
equation (Vt(0) = 0) and shape parameters or dis- 
tortions of gravitationally lensed images (VVr(0)), 
all at values corresponding to the observed image 
locations, are considered as rigid linear constraints on 
the mass distribution. 

Writing the lens equation at an observed image lo- 
cation 0\ 



(3 = 0, 



Dds 



^2 K mn^0lpmn{Ol), (2-5) 



each image supplies us with two-component con- 
straint equation. But we have to solve for the un- 
known source position j3. Thus the number of con- 
straints on the mass distribution from multiple im- 
ages is 2 ((images) — (sources)). These strong lensing 
constraints are linear equality constraints. 




Origin 

Figure 1 . Illustration of the notation used in the text for 
explaining the constraints from weak lensing. The X >, 9 y i 
axes are the principal axes of the observed stretched ar- 
clet, and are rotated by an observed angle <p with respect 
to the main coordinate axes 9 x ,0 y . 



The constraints from the observed shape param- 
eters of distorted images (arclets) can be obtained 
from the the Hessian of the arrival time function. If 
an arclet at 0\ is observed to be stretched along the 
direction 6 y / by at least a factor of e compared to 
the perpendicular direction 6 x i (see Fig. 1), then by 
considering the 9 x i, 9 y i components of the inverse am- 
plification matrix we have 
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If we can independently infer the image parity we can 
remove the absolute value signs in equation ( |2-6| ) , 
which will then supply us with a linear constraint 
equation on the because rotation to 9 x i,9 y i co- 
ordinates is a linear transformation. 

While constraints of the form ( |2-6| ) can be used for 
arclets in multiply-imaged systems (and we do so), 
they are most useful for singly imaged arclets (weak 
lensing). For statistical distortions the form of the 
constraint stays exactly the same. If the shear field is 
known accurately enough the inequality can become 
an equality, and can be supplemented by the addi- 
tional constraint 
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In our work the weak lensing constraints are ex- 
pressed in implicit form (Eq. 2-6), whereas in previ- 
ous work (e.g. KS93, Kaiser 1995, Seitz & Schneider 
1996) it is usual to express these in explicit form. The 
advantage of the implicit form is that the constraints 
remain linear in all lensing regimes. 

2.3. Breaking the mass sheet degeneracy 

The mass sheet degeneracy is as follows. Consider 
the lens equation in the form 



(3 = 
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ds 
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simply amounts to multiplying both sides of Eq. (2 
by r. Clearly, image positions will be unaffected by 
such a transformation. Also, since the magnifica- 
tion matrix is just 86 / d(3 the transformation will just 
rescale the whole matrix by r _1 ; in other words, rel- 
ative shear and magnifications will be unaffected but 
absolute amplification will change. Physically, the 
transformation (2-9) amounts to rescaling the mass 
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everywhere and then adding a constant mass sheet; 
hence the name mass-sheet degeneracy. 

It is sometimes thought that measuring absolute 
amplifications is the only way to break the mass-sheet 
degeneracy. This is not so. It is easy to verify that 
the transformation ( |2-9| ) will not have the same in- 
variance property for two different source redshifts at 
the same time. Thus lensing data with more than one 
source redshift breaks the mass-sheet degeneracy. 



2.4. Mass Profile Reconstruction 

The above lensing constraints, together with the 
physical requirement of n mn > 0, by themselves are 
insufficient to constrain the lens uniquely. The fact 
that the number of lensing constraints is far less than 
the number of pixels leaves a vast family of mass dis- 
tributions, all perfectly consistent with the lensing 
observations. It is now necessary to add more infor- 
mation, based on some criteria for physical plausibil- 
ity. 

As discussed in Paper I, it is particularly interesting 
to consider mass distributions that minimize 
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where L mn represents the light profile normalized to 
unit total luminosity. In Eq. fl2-10|) , the first term 
tends to minimize the M/L variations while the sec- 
ond term smoothes the mass map on scales of < a. 
Such minimization may be considered as regulariza- 
tion with respect to the light distribution and the 
smoothing scale a, and is readily implemented via 
standard numerical algorithms, such as the NAG rou- 
tine E04NFF. It serves two purposes: (i) as a test 
whether light is indeed a fair or biased tracer of mass, 
and (ii) a basis for Monte-Carlo simulations to esti- 
mate uncertainties in the mass map. We go into de- 
tails in the next section, but for now we emphasize 
strongly that our technique considers the light distri- 
bution as a secondary information subservient to the 
rigid constraints from the observed lensing data. 



3. LENSED OBSERVABLES IN ABELL 2218 

Abell 2218 is an exceptionally rich lensing cluster, 
at redshift z = 0.175, which hosts 7 multiple image 
systems and over 100 arclets. In the present paper, 
we will use all the secure informations that lensing 
provides in this first application of combined strong 
and weak lensing. Our analysis of A2218 is based on 
the archival WFPC2/HST images and we take the 
redshifts of the resolved images from Ebbels et al. 
(1998). 



Throughout this paper we use Q = 1; A = 0; Hq = 
50h km/s/Mpc which corresponds to an angular size 
of 1" = 3.84 kpc at the redshift of Abell 2218. 

To refer to individual cluster and background 
galaxies, we will follow the three-digit numbering 
scheme of Le Borgne et al. (1992); but for objects 
referred to often we adopt simpler names below. 

3.1. Multiply Lensed features 

Of the multiply imaged background galaxies, we 
use the three that have secure spectroscopic redshifts 
and one with a photometric redshift. Positions and 
redshifts are listed in Table 1 and plotted in Fig. 2. 



ID 


z 


Image 


X 


V 


384 


2.515 


al 


14.5 


19.1 




2.515 


a2 


17.0 


15.4 


328 


0.702 


bl 


-14.7 


16.1 


389 


0.702 


b2 


-17.8 


-15.4 


289 


1.034 


cl 


-61.7 


7.2 




1.034 


c2 


-61.5 


0.6 


730 


1.1 


dl 


-75.2 


-1.4 




1.1 


d2 


-74.6 


-6.0 




1.1 


d3 


-73.6 


-9.0 



Table 1: Positions of various sets of multiple images (x,y 
in arcsec), as taken from the HST image, used in our re- 
construction. 
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Figure 2. The positions of the multiple images and the 
arclets in the WFPC2/HST. 



• The 'a' system (z = 2.515j: 
The spectacular arc system #384 clearly shows an in- 
ternally symmetric pattern of unresolved knots which 
characterizes the system as a fold arc comprising two 
merging images (which we will call al & a2) with 
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opposite parities. It is identified to be an image of 
a background star forming galaxy, with the visible 
knots representing the H II regions. 

• The 'b' system (z = 0.702j: 

The identification of the multiple images in this sys- 
tem is somewhat complicated. The red arc #359 was 
initially identified by Pello et al. (1992) to be an im- 
age of a background spheroidal galaxy. It was later 
interpreted by Kneib et al. (1995) using a ground- 
based image to be two merging images forming a fold 
arc with #328 as a counter image. However, the im- 
proved resolution of HST revealed that such a config- 
uration for the arc #359 cannot be true (KESCS96) 
since a faint extension connecting it to #337 is now 
revealed, thus strongly suggesting it as a counter im- 
age too. On the other hand, #337 and #389 yield 
a similar color to #359, thus indicating that #389 is 
another counter image of #359. Since such a system 
is complicated and no simple model can explain it, 
we tried various input configurations. The only plau- 
sible solutions can be obtained when #328 and #389 
(which we call bl and b2) are considered as two im- 
ages from the same source, and #359 (which we call 
'b3?') is an arclet arising from another source, possi- 
bly a different component of the same lensed galaxy. 

• The 'c' system (z = 1.034): 

The blue arc #289 has a large amount of internal 
structure in the southern part, which is highly mag- 
nified but appears to be slightly distorted and thus 
characterized to be singly imaged while the tail visible 
on the north is is highly elongated with tiny breaks 
and definitely multiply imaged. Nevertheless, this 
complex configuration can be explained by a back- 
ground source lying near a cusp caustic but with a 
high portion of it being outside the cusp and ob- 
served as the southern end of the arc, while the por- 
tion within the caustic is observed as the northern 
end. We consider the northern part as two merging 
segments (and call these cl and cl), and consider the 
southern part (which we call c3?) as a singly imaged 
arclet. 

• The 'd' system (z = 1.1 ± 0.3J: 

This is a faint thin arc, known as #730, with a num- 
ber of bright knots visible. Three elongated compo- 
nents (which we call dl, d2, and d3) can be easily 
distinguished along the arc which characterizes it as 
a cusp arc. 

3.2. Secure single distorted arclets: 

Ebbels et al. (1998) have spectroscopically mea- 
sured the redshifts of a number of arclets with differ- 
ent degrees of accuracy in the HST/WFPC2 area. We 



use all the arclets for which they have quality- 1 red- 
shifts, and arclets with quality-2 redshifts in regions 
where quality- 1 redshifts are not available. Table 2 
gives details of the 18 arclets we use. The last two 
arclets constitute parts of multiple image systems. 
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y 


6 


e 


145 


0.628 


—2.5 


103.9 


26.57 


2.5 


159 


0.564 


5.1 


105.4 


80.54 


2.0 


158 


0.723 


—41.8 


62.9 


19.65 


2.2 


242 


0.635 


30.3 


88.8 


117.9 


5.0 


231 


0.563 


-9.7 


57.8 


29.48 


4.0 


381 


0.521 


35.9 


34.6 


157.4 


2.5 


317 


0.474 


-7.3 


27.5 


167.7 


1.5 


306 


0.450 


-97.5 


-44.8 


135.0 


1.8 


464 


0.476 


10.6 


-23.3 


45.0 


6.0 


467 


0.475 


17.0 


-19.7 


48.01 


4.0 


297 


0.450 


-21.7 


23.2 


65.22 


5.0 


456 


0.538 


20.6 


-11.6 


49.76 


3.5 


238 


0.635 


-44.4 


20.1 


49.40 


4.0 


273 


0.800 


-41.4 


16.6 


68.75 


2.5 


205 


0.693 


-72.1 


16.3 


90.00 


1.2 


431 


0.675 


39.9 


21.9 


110.0 


6.0 


444 


1.030 


24.6 


-2.8 


75.96 


7.0 


359 (b3?) 


0.702 


-20.5 


-3.7 


98.75 


5.0 


289 (c3?) 


1.034 


-63.7 


-8.9 


77.20 


6.0 



Tabic 2: Observational parameters of the arclets used 
in our reconstruction. The positions x and y are in arcsec 
and the position angle 6 is in degrees, while e is a lower 
bound on the elongation. The last two arclets, b3? and 
c3? are probably parts of the corresponding multiple im- 
age systems in Table 1, but we have not used these as 
multiple-image constraints. We have taken all the arclets 
as minima of the arrival time, except for b3? which is 
evidently a saddle. 



3.3. Luminosity Distribution 

The apparent magnitudes in the i?-band of the 
foreground galaxies are taken from Le Borgne et al. 
(1992). We considered only those with apparent mag- 
nitudes < 20. and we find that there are 33 of them 
enclosed within the HST image. The pixellated light 
distribution L mn is obtained by replacing each of the 
33 galaxies by a Gaussian light profile of dispersion 
10". Fig. 3 shows a contour plot of this smoothed 
luminosity distribution. 
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Figure 3. The luminosity distribution of the 

HST/WFPC2 region of A2218. The positions of the 33 
galaxies used are marked with filled circle. 



Figure 4. Mass distribution of Abell 2218, reconstructed 
using strong and weak lensing constraints, and regulariz- 
ing with respect to the light distribution and a smoothing 
scale of 6". The density is in units of the critical density 
for sources at infinity, £ C rit = 3.1 x 1O 14 /i 5O 1 M0 arcsec -2 ; 
contours are in steps of 0.1. Later figures in this paper 
always use the same units, and use the same contour steps 
unless noted otherwise. The total mass in the field is 
2.75 x W 14 h^M & . 



It is of interest to examine also other reconstructed 
mass maps obtained using different criteria. For ex- 
ample, one can reconstruct a mass map by assum- 
ing zero light distribution for the cluster (i.e., setting 
Lmn = in Eq. 2-10| ); this is shown in Fig. 5, and 
is roughly the minimal matter distribution needed to 
reproduce the data. Or one can reconstruct a mass 
map with a much higher smoothing parameter; Fig. 6 
shows a mass map using a = 19". It is reassuring to 
see that all the main features are still present in these 
cases. 



4. RECONSTRUCTED MASS DISTRIBUTIONS AND 
UNCERTAINTIES 

Our optimal mass map of A2218 is shown in Fig. 4. 
It uses a pixel size of 2.29" and regularizes with re- 
spect to the light distribution with smoothing scale 
g = 6", while strictly obeying all the strong and 
weak lensing constraints. (We found a = 6" to be 
the smallest value that eliminated obvious artifacts 
in the mass map. We use this value throughout this 
paper except where noted otherwise.) Our mass map 
is broadly similar to the previous parametric mod- 
els (K95, KESCS96), which show that the cluster is 
clearly bimodal with one clump being less massive, 
and with detailed substructures implied by the lens- 
ing data. 
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Figure 5. Mass distribution reconstructed for the zero- 
light case. The total mass in this case is 1.06x 10 14 ft-^ 1 M Q . 




Figure 6. Mass map using much more smoothing: 
a = 19". 

The spikes visible in all our mass maps are due to 
fact we are using local constraints from sparsely sam- 
pled data, i.e., arclets. The peaks of the spikes are 
very robust between different reconstructions, but the 
wings are very variable. As a result the total mass 
shows large variation, though the allowed variation is 
quite well bounded above and below. Statistical dis- 
tortion maps going to larger radii would, we expect, 
constrain the total mass much better and hence re- 
duce this problem. We leave this extension for future 
work, hopefully with several HST pointings. 
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Figure 7. Contour map of the projected mass distribu- 
tion. This corresponds to Fig. 4 and is plotted here for 
comparison with Fig. 8. 
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Figure 8. Fractional uncertainty of the mass distribu- 
tion. The crosses mark the position of the two dominant 
cD galaxies in A2218. 

Our reconstruction technique naturally lends itself 
towards calculating error estimates on the pixellated 
mass distribution. We use a Monte-Carlo procedure, 
randomizing L mn by re-shuffling the positions of the 
cluster galaxies or by rotating the entire light distri- 
bution by various angles, and constructing an ensem- 
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ble of new mass maps. The fractional dispersion 

( K mn) 



( K r, 



1 



(4-1) 



over this ensemble is a measure of the pixel-by-pixel 
fractional uncertainty in the mass map. This is shown 
in Fig. 8 as a contour plot of the computed An mn . 
Regions in the lens plane enclosed by A« mn < 0.2, are 
very well constrained. Clearly, the best constrained 
regions have a high number density of images. In 
Fig. 9 we quantify this by plotting AK mn binned over 
circles of radius 10" versus the the number of images 
enclosed. As expected, the fractional uncertainty de- 
creases with the number of images. 
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Figure 9. The fractional uncertainty versus the num- 
ber of images. Each dot represents a circle of radius 10" 
over which AK mn and the number of images have been 
averaged. The dashed curve is a linear fit. 



5. MASS-LIGHT DISCREPANCIES 

While our mass reconstruction broadly resembles 
previous parametric models, we find significant differ- 
ences that contradict a basic assumption of the para- 
metric modelling method, that of modelling the clus- 
ter as a smooth mass distribution plus small clumps 
associated with bright galaxies. Two of these differ- 
ences are worth discussing in detail. 

First, we find that the highest mass peak is signifi- 
cantly offset from the brightest light peak, the center 
of the cD galaxy #391. The offset is highly signifi- 
cant given the uncertainties. To test this further, we 
did a mass reconstruction in a smaller field around 
this main peak but with a finer pixel size of 1". We 
find (see Fig. 10) that the offset is still persistent and 



is about 10" towards the second peak of the cluster. 
It is interesting to note that the direction of the off- 
set between the mass and light peaks detected from 
our modelling is similar to that found by Markevitch 
(1997) between the temperature and light peak. 




-20 \ 



X[arcsec] 



Figure 10. Mass reconstruction of a smaller field, cen- 
tered on the cD galaxy #391, with 1" pixels; a cross marks 
the cD center. 



-20 




X[arcsec] 

Figure 11. Mass reconstruction of a smaller field cen- 
tered on the galaxy #244, with 1" pixels; crosses mark 
this galaxy, and two other fainter galaxies. 



Second, we find that the secondary mass peak is 
not centered on the bright galaxy #244. Again, we 
did a mass reconstruction on a smaller field using 1" 
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pixels, this time near #244; this is shown in Fig. 11. 
The mass distribution in this region, has two peaks in 
the directions of two galaxies at least 2.5 mag fainter 
in R than #244, while #244 itself is in a valley be- 
tween these two peaks. This is also highly significant 
given the uncertainties, and is seen both in the main 
reconstruction and in the 'zoom-in'. If the three mass 
concentrations are indeed associated with the identi- 
fied galaxies, the 2.5 mag difference implies that the 
associated M/L ratios vary by more than a factor of 
10. 

These new results we reported are robust to pixel 
size, shape and to the light distribution of the cluster. 

6. COMPARISON WITH X-RAY MODELS 

The X-ray emission by clusters is usually modelled 
assuming that the underlying potential of the cluster 
is spherically symmetric and that the hot gas is in 
hydrostatic equilibrium. It is important to test these 
assumptions by comparing with lensing, which does 
not depend on these assumptions. Miralda-Escude 
& Babul (1995) argued that the strong lensing data 
could not be reconciled with an equilibrium spherical 
model for the the X-ray emitting has at the observed 
temperature — to agree with lensing the gas temper- 
ature would have to be much higher, or equivalently, 
lensing required at least twice as much mass as the 
gas model. They suggested that the gas might be 
partly supported by turbulent motion or by magnetic 
fields, or might be multi-phased. Another possibility 
is that the cluster is still undergoing a merging phase, 
and hence the gas would not be expected to trace the 
gravitational potential. 
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Figure 12. Projected mass within different radii from the 
main cD galaxy #391, from an X-ray model by Allen et 
al. 1998 and from lensing. The two lensing reconstruc- 
tions indicated correspond to Figs 4 (maximum) and 5 



(minimum) . 

In Fig. 12 we plot the enclosed projected mass from 
a model for the hot gas by Allen (1998) and enclosed 
mass from our reconstructions. The lowest enclosed 
mass from lensing we got for the 'zero-light' recon- 
struction shown in Fig. 5; the highest is for our opti- 
mal reconstruction (Fig. 4) which tends to minimize 
M/L variation and thus extrapolates mass into re- 
gions with now lensing data; other models correspond 
to various other regularizations. We see from Fig. 12 
that the mass discrepancy is a factor of 2.5 for even 
the unrealistic-looking model of Fig. 5. 

7. ANALYSIS OF INDIVIDUAL MULTIPLE IMAGE 
SYSTEMS 

We have examined the arrival time surfaces, crit- 
ical curves, and caustics for all the multiple image 
systems and some of the arclets. Arrival time con- 
tours are a good way of verifying a mass distribution 
reproduces the given image position properties. In 
this section we discuss two of the multiple image sys- 
tems in detail, and report predicted counter- images, 
if any. 



7.1. The al-a2 system 




X[orcsec] 

Figure 13. Arrival time surface for the al-a2 system at 
z=2.515. Filled circles mark the observed positions of al 
and a2. A third image is predicted roughly at (—5", —23"). 

The arrival time surface (see Fig. 13) reproduces 
the image positions of the two segments al and a2 
comprising the fold arc. Moreover, the arrival time 
surface predicts a third image 20" below the domi- 
nant mass clump, very close to the observed arclet 
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#468, also predicted by K96 counter image for 
the same arc. The mirror symmetry seen across the 
arc al-a2 implies that the critical curve is passing 
through them. Exploring the critical curves at the 
redshift (z = 2.515) of the system confirms that al- 
a2 is indeed a fold arc resulting from a background 
source lying close to a beak-to-beak caustic. The crit- 
ical curves and caustics implied by our mass map are 
shown in Fig. 14. 
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Figure 14. The critical curves and caustics at the redshift 
of the al-a2 system (z = 2.515). Fhc upper panel shows 
the critical curves with the two filled circles marking the 
positions of the two merging images. The lower panel is 
the corresponding caustics with the filled circle marking 
the predicted source position. 



7.2. The bl-b2 system 

The arrival time surface (see Fig. 15) in addition to 
reproducing the positions of bl and b2 as required, 
also predicts an image near the observed position of 
b3?, but not exactly coinciding with it, which sug- 
gests that b3? is an image of a different component 
of the same lens galaxy. Fig. 16 shows the critical 
curves caustics at the redshift (z = 0.702) of this sys- 
tem. The predicted position of the source is in a re- 
gion where the caustics are very convoluted and this 
would result in a complex image configuration. As an 
auxiliary test, we inspect the arrival time surface for 
the arclet b3? itself (see Fig. 17). The position of the 
arclet b3? is of course exactly reproduced and, more- 
over, four extra images are predicted. The position of 
two of the predicted images matches fairly with those 
of the bl, and b2. This result strongly supports the 
scenario that the images bl, b2 and #337 are counter 



images of the arclet b3?. 
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Figure 15. Arrival time surface for the bl-b2 system at 
z=0.702. Filled circles mark the observed positions of bl 
and b2. 
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Figure 16. The critical curves and caustics at redshift 
z = 0.702. Fhe critical curve (upper panel) with filled cir- 
cles marking positions of observed images. Caustics (lower 
panel) with filled circle marking the predicted position of 
the source. 
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Figure 17. Arrival time surface for b3? at z=0.702. The 
filled circle marks the observed position of b3?. 



7.3. Statistical Magnification Map 

Since most work on cluster lensing at present is 
based on statistical distortion of images, it is interest- 
ing to examine the statistical magnification map. We 
calculate the statistical magnification as follows. We 
first divide the field into square regions of side 10" or 
20". Within each square we compute the magnifica- 
tion matrix at 25 random points, using our optimal 
mass reconstruction and redshifts randomly chosen 
from those of the multiply imaged systems; of these 
we discard any points corresponding to axis ratios 
> 5 and average the rest. This averaged matrix we 
call the statistical magnification. Such a procedure 
roughly mimics the observational procedure of av- 
eraging ellipticities while discarding obvious arclets, 
though of course it is not the same because for real 
data the absolute magnification is not usually avail- 
able. 

Fig. 18 shows the statistical magnification we ob- 
tained by averaging over 10" and 20" squares. It illus- 
trates a cautionary fact: in the strong lensing region 
the statistical magnification is dominated by noise. 
The reason is that a critical curve might be crossed 
and this can change the shear from ~ 1 to oo and 
back to ~ 1 again over a 10" scale. Discarding arclets 
with high axis ratios does not cure this problem — the 
left panel in Fig. 18 shows two places with statistical 
magnification corresponding to axis ratios > 5 even 
though individual location with such high axis ratios 
have been excluded; this is because two low-axis-ratio 
magnification matrices on opposite sides of a critical 
curve can result in a high-axis-ratio one if averaged. 
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Figure 18. Statistical magnification maps from our opti- 
mal mass reconstruction, averaging over 10" squares (top), 
and 20" squares (bottom). The size of each ellipse repre- 
sents the convergence, the shape represents the shear. 

We conclude that statistical shear must not be used 
in the strong lensing region of clusters; instead, con- 
straints from each image should be used separately. 

8. CONCLUSIONS 

The present paper is, to our knowledge, the first to 
combine constraints from strong and weak lensing in 
cluster mass reconstruction. The method we describe 
recovers a pixellated mass distribution that strictly 
obeys constraints from lensing observations and uses 
the light distribution of the cluster as subordinate 
information that may be overridden by lensing con- 
straints. We explore the projected (dark) matter dis- 
tribution of Abell 2218 with an unprecedented level 
of detail. Our mass map shows that the primary mass 
peak is offset from the the light peak and that pro- 
jected mass-to-light variations of galaxy-sized compo- 
nents are severely inconsistent with the galaxy M/L 
scaling deduced from dynamical arguments. In gen- 
eral, these imply that mass does not follow light on 
a range of length scales in A2218. We also confirm 
and elaborate a previous result that current X-ray 
mass models significantly underestimate the mass of a 
cluster that is able to reproduce the observed lensing. 
Our results conclude that mass estimates from lensing 
are at least 2.5 times that from X-ray models. This 
suggests that in at least some clusters intra-cluster 
hot gas does not trace the gravitational potential. 
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APPENDIX 
INTEGRALS OVER GAUSSIAN PIXELS 



Integrating (2-4) yields the coefficient of the mn-th pixel's contribution to the potential: 

a 2 26 2 
ipmn(0x, Oy) = — [In (u) - Ei(-u) +je], u = 



(Al) 



where Ei denotes the exponential integral and 7^ is Euler's constant. Thus the first derivatives of the coefficients 
of the deflection potential, i.e., the components of the deflection angle, are 



[1 - exp(-u)] 
[1 - exp(-u)] . 



(A2) 



The second derivatives of the potential are 

dxx1pmn{@xi$y) — 

dyy^mni^xiGy) — 
dxy1^mn{@xi@y) — 
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[1 - exp(-n)] + exp(-n), 
26 2 

[1 - exp(-n)] + exp(-n), 



1 



exp(— u) [1 — exp(— u)] 
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(A3) 



